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Abstract 


A  fully  analytical  model  for  the  prediction  of  the  one-dimensional  mean  wind  speed,  kine¬ 
matic  shear  stress,  turbulence  kinetic  energy,  and  velocity  variances  in  a  horizontally  homo¬ 
geneous  canopy  is  described.  The  basis  of  the  model  centers  around  an  analytical  solution 
obtained  by  Massman  and  Weil  (Boundary-Layer  Meteorology,  91,  pp.  81-107,  1999)  for  the 
turbulence  kinetic  energy  derived  from  a  particular  form  of  a  second-order  closure  model  for 
canopy  flows.  Using  this  analytical  model  for  the  turbulence  kinetic  energy,  it  is  shown  how 
to  derive  analytical  expressions  for  the  within  canopy  velocity  variances  using  two  meth¬ 
ods:  namely,  (1)  equilibrium  partitioning  whereby  the  velocity  variances  are  assumed  to  be 
proportional  to  the  turbulence  kinetic  energy,  with  the  proportionality  constants  the  same 
as  those  at  the  canopy  top;  and,  (2)  algebraic  stress  model  partitioning  derived  from  an 
algebraic  stress  model  version  of  an  existing  second-order  closure  model  for  flows  proposed 
by  Launder,  Reece,  and  Rodi  (Journal  of  Fluid  Mechanics,  68,  pp.  537-566,  1975)  that  has 
been  modified  to  include  the  effects  on  the  turbulence  of  the  extra  physics  arising  from  the 
interaction  of  the  airflow  with  the  obstacles  in  the  array  (viz.,  the  form  drag  on  the  canopy 
elements,  and  the  role  of  the  momentum  absorption  over  an  extended  volume  of  space  on 
the  turbulent  motions  within  the  canopy). 

The  predictions  from  the  analytical  model  have  been  compared  with  some  measurements  of 
mean  flow  and  turbulence  statistics  obtained  in  two  urban  canopies:  namely,  (1)  Tombstone 
canopy;  and,  (2)  arrays  of  regularly  spaced  cubical  obstacles  at  various  plan  area  indices. 
It  is  shown  that  the  model  is  moderately  successful  at  providing  estimates  of  a  number 
of  important  flow  quantities;  namely,  mean  wind  speed  and  various  turbulence  statistics. 
Given  the  limited  data  that  was  available,  it  has  not  been  possible  to  ascertain  whether  the 
algebraic  stress  model  partitioning  of  the  turbulence  kinetic  energy  offers  any  significant 
advantages  in  the  determination  of  the  component  velocity  variances  over  the  simpler  equi¬ 
librium  partitioning.  With  the  testing  of  the  analytical  model  against  the  urban  canopy 
flows  and  its  modest  success,  it  is  suggested  that  this  analytical  model  for  canopy  flows  can 
probably  be  applied  unmodified  to  a  wide  range  of  canopies,  both  vegetative  and  urban. 
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Executive  summary 


Introduction:  Incidents  in  Tokyo  involving  the  use  of  chemical  and  biological  warfare 
(CBW)  agents  by  a  terrorist  organization  have  mobilized  both  military  and  public  concerns 
on  the  use  of  these  weapons  against  military  and  civilian  populations  in  built-up  or  urban 
areas  consisting  of  large  groups  of  obstacles.  In  response  to  these  increasing  concerns,  the 
prediction  of  the  dispersion  of  CBW  agents  released  within  an  urban  area  has  recently  be¬ 
come  the  subject  for  many  investigations  since  information  obtained  from  urban  dispersion 
models  will  be  used  to  form  the  basis  on  which  government  and  military  commanders  (as 
well  as  other  decision  makers)  will  respond  and  act  to  mitigate  the  effects  of  a  particular 
incident.  Physically-based  models  such  a  Lagrangian  Stochastic  (LS)  models  exist  for  the 
prediction  of  plume  or  cloud  dispersion  in  complex  environments,  and  certainly  can  be  ap¬ 
plied  to  the  problem  of  urban  dispersion.  The  prediction  of  dispersion  using  LS  models, 
which  involves  the  calculation  of  the  trajectories  of  a  large  number  of  tagged  fluid  particles 
representing  the  contaminant  mass  and  the  derivation  of  the  ensemble-average  concentra¬ 
tion  from  the  calculated  particle  distribution,  is  enticingly  simple  in  concept  and  general  in 
approach.  However,  for  urban  dispersion,  LS  models  require  to  be  provided  information  on 
the  mean  and  statistical  properties  of  the  turbulent  wind  field  within  and  above  the  urban 
canopy  (e.g.,  mean  wind  speed,  shear  stress,  velocity  variances,  etc.). 

The  mean  flow  and  turbulence  within  and  immediately  above  an  urban  canopy  (wherein  the 
atmospheric  flow  is  influenced  by  the  individual  roughness  elements  or  obstacles)  possess  a 
fully  three-dimensional  structure.  Although  the  development  of  three-dimensional  dynami¬ 
cal  models  of  the  mean  flow  and  turbulence  in  simplified  urban  canopies  consisting  of  large 
groups  of  regularly  spaced  obstacles  is  progressing,  these  models  tend  to  require  extensive 
computational  resources.  These  flow  models  involve  the  solution  of  the  conservation  equa¬ 
tions  of  mass,  momentum,  and  energy,  and  the  detailed  information  obtained  from  these 
models  can  be  used  to  “drive”  physically-based  dispersion  models  to  provide  estimates  of 
the  contaminant  concentrations.  However,  currently  the  computational  demands  of  these 
three-dimensional  models  will  preclude  their  use  for  emergency  response  situations  which 
require  the  ability  to  generate  a  dispersion  prediction  in  a  time  frame  that  will  permit  pro¬ 
tective  actions  to  be  taken  (i.e.,  time  scales  on  the  order  of  seconds  to  minutes).  In  view 
of  this,  we  argue  that  for  many  practical  applications  it  is  convenient  to  consider  spatial 
averages  of  the  mean  wind  and  turbulence  within  an  urban  canopy  obtained  by  averaging 
horizonally  the  mean  wind  and  turbulence  statistics  over  an  area  that  is  much  greater  than 
the  spacings  between  the  individual  roughness  elements  comprising  the  urban  canopy.  This 
provides  a  simplified  one-dimensional  mean  wind  and  turbulence  model  for  urban  canopy 
flow  whereby  the  flow  statistics  depend  only  on  the  height  z.  In  this  simplified  description 
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of  urban  flow,  the  spatially-averaged  mean  wind  and  turbulence  vary  only  with  height,  and 
the  detailed  three-dimensional  description  of  flow  and  turbulence  characteristics  about  the 
individual  buildings  is  lost. 

The  primary  objective  of  this  report  is  to  describe  and  test  a  simple,  analytical  one¬ 
dimensional  model  for  the  prediction  of  mean  flow  and  turbulence  within  and  above  an 
urban  canopy.  This  model  was  developed  originally  for  vegetative  canopies  (e.g.,  forests, 
cereal  crops,  etc.),  but  will  be  applied  here  to  urban  canopies.  While  there  are  similarities 
between  vegetative  and  urban  canopies,  there  are  enough  important  differences  between 
these  two  types  of  canopies  (e.g.,  characteristic  density  of  roughness  elements  and  their 
stiffness  or  ability  to  react  to  the  drag  forces)  to  require  an  independent  evaluation  of 
the  model  for  the  urban  canopy  application.  It  is  hoped  that  this  simple,  analytical  one¬ 
dimensional  model  for  urban  canopy  flows  can  be  used  to  provide  useful  input  information 
on  the  statistics  of  the  wind  field  required  to  drive  LS  models  of  urban  dispersion. 

Results:  A  fully  analytical  model  for  the  prediction  of  the  one-dimensional  mean  wind 
speed,  kinematic  shear  stress,  turbulence  kinetic  energy,  and  velocity  variances  in  a  hori¬ 
zontally  homogeneous  canopy  is  described.  The  basis  of  the  model  centers  around  an  analyt¬ 
ical  solution  obtained  by  Massman  and  Weil  (Boundary-Layer  Meteorology,  91,  pp.  81-107, 
1999)  for  the  turbulence  kinetic  energy  derived  from  a  particular  form  of  a  second-order  clo¬ 
sure  model  for  canopy  flows.  Using  this  analytical  model  for  the  turbulence  kinetic  energy, 
it  is  shown  how  to  derive  analytical  expressions  for  the  within  canopy  velocity  variances 
using  two  methods:  namely,  (1)  equilibrium  partitioning  whereby  the  velocity  variances 
are  assumed  to  be  proportional  to  the  turbulence  kinetic  energy,  with  the  proportionality 
constants  the  same  as  those  at  the  canopy  top;  and,  (2)  algebraic  stress  model  partitioning 
derived  from  an  algebraic  stress  model  version  of  an  existing  second-order  closure  model  for 
flows  proposed  by  Launder,  Reece,  and  Rodi  (Journal  of  Fluid  Mechanics,  68,  pp.  537-566, 
1975)  that  has  been  modified  to  include  the  effects  on  the  turbulence  of  the  extra  physics 
arising  from  the  interaction  of  the  airflow  with  the  obstacles  in  the  array  (viz.,  the  form 
drag  on  the  canopy  elements,  and  the  role  of  the  momentum  absorption  over  an  extended 
volume  of  space  on  the  turbulent  motions  within  the  canopy). 

The  predictions  from  the  analytical  model  have  been  compared  with  some  measurements  of 
mean  flow  and  turbulence  statistics  obtained  in  two  urban  canopies:  namely,  (1)  Tombstone 
canopy;  and,  (2)  arrays  of  regularly  spaced  cubical  obstacles  at  various  plan  area  indices. 
It  is  shown  that  the  model  is  moderately  successful  at  providing  estimates  of  a  number 
of  important  flow  quantities;  namely,  mean  wind  speed  and  various  turbulence  statistics. 
Given  the  limited  data  that  was  available,  it  has  not  been  possible  to  ascertain  whether  the 
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algebraic  stress  model  partitioning  of  the  turbulence  kinetic  energy  offers  any  significant 
advantages  in  the  determination  of  the  component  velocity  variances  over  the  simpler  equi¬ 
librium  partitioning.  With  the  testing  of  the  analytical  model  against  the  urban  canopy 
flows  and  its  modest  success,  it  is  suggested  that  this  analytical  model  for  canopy  flows  can 
probably  be  applied  unmodified  to  a  wide  range  of  canopies,  both  vegetative  and  urban. 


Significance  and  Future  Plans:  The  one-dimensional  analytical  model  for  canopy  flows 
considered  in  this  report  provides  all  the  statistics  of  the  wind  field  (mean  wind  speed, 
kinematic  shear  stress,  velocity  variances,  and  viscous  dissipation  rate)  required  as  input  to 
drive  a  physically-based  model  for  the  prediction  of  urban  dispersion.  These  flow  statistics 
can  be  calculated  rapidly  using  the  analytical  model  in  comparison  to  the  more  compu¬ 
tationally  intensive  second-order  closure  models  for  canopy  flows.  Taken  as  a  whole,  this 
simple  canopy  flow  model  can  probably  be  used  to  “drive”  a  Lagrangian  stochastic  disper¬ 
sion  model  for  the  prediction  of  dispersion  of  materials  released  within  a  canopy.  The  next 
logical  step  in  this  work  would  be  to  construct  a  Lagrangian  stochastic  model  satisfying  the 
well-mixed  criterion  that  uses  the  simple  analytical  model  to  provide  the  input  mean  flow 
and  turbulence  quantities  required. 


Yee,  E.,  2000.  A  One-Dimensional  Mean  Wind  and  Turbulence  Model  for  a  Uniform  Urban 
Canopy.  DRES  TR  2000-071,  Defence  Research  Establishment  Suffield. 
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Introduction 


Incidents  in  Tokyo  involving  the  use  of  chemical  and  biological  warfare  (CBW)  agents  by 
a  terrorist  organization  have  mobilized  both  military  and  public  concerns  on  the  use  of 
these  weapons  against  military  and  civilian  populations  in  built-up  or  urban  areas  consist¬ 
ing  of  large  groups  of  obstacles.  In  response  to  these  increasing  concerns,  the  prediction 
of  the  dispersion  of  CBW  agents  released  within  an  urban  area  has  recently  become  the 
subject  for  many  investigations  since  information  obtained  from  urban  dispersion  models 
will  be  used  to  form  the  basis  on  which  government  and  military  commanders  (as  well  as 
other  decision  makers)  will  respond  and  act  to  mitigate  the  effects  of  a  particular  incident. 
Physically-based  models  such  a  Lagrangian  Stochastic  (LS)  models  exist  for  the  prediction 
of  plume  or  cloud  dispersion  in  complex  environments,  and  certainly  can  be  applied  to  the 
problem  of  urban  dispersion.  The  prediction  of  dispersion  using  LS  models,  which  involves 
the  calculation  of  the  trajectories  of  a  large  number  of  tagged  fluid  particles  representing 
the  contaminant  mass  and  the  derivation  of  the  ensemble-average  concentration  from  the 
calculated  particle  distribution,  is  enticingly  simple  in  concept  and  general  in  approach. 
However,  for  urban  dispersion,  LS  models  require  to  be  provided  information  on  the  mean 
and  statistical  properties  of  the  turbulent  wind  field  within  and  above  the  urban  canopy 
(e.g.,  mean  wind  speed,  shear  stress,  velocity  variances,  etc.). 

The  mean  flow  and  turbulence  within  and  immediately  above  an  urban  canopy  (wherein  the 
atmospheric  flow  is  influenced  by  the  individual  roughness  elements  or  obstacles)  possess  a 
fully  three-dimensional  structure.  Although  the  development  of  three-dimensional  dynami¬ 
cal  models  of  the  mean  flow  and  turbulence  in  simplified  urban  canopies  consisting  of  large 
groups  of  regularly  spaced  obstacles  is  progressing,  these  models  tend  to  require  extensive 
computational  resources.  These  flow  models  involve  the  solution  of  the  conservation  equa¬ 
tions  of  mass,  momentum,  and  energy,  and  the  detailed  information  obtained  from  these 
models  can  be  used  to  “drive”  physically-based  dispersion  models  to  provide  estimates  of 
the  contaminant  concentrations.  However,  currently  the  computational  demands  of  these 
three-dimensional  models  will  preclude  their  use  for  emergency  response  situations  which 
require  the  ability  to  generate  a  dispersion  prediction  in  a  time  frame  that  will  permit  pro¬ 
tective  actions  to  be  taken.  In  view  of  this,  we  argue  that  for  many  practical  applications 
it  is  convenient  to  consider  spatial  averages  of  the  mean  wind  and  turbulence  within  an 
urban  canopy  obtained  by  averaging  horizonally  the  mean  wind  and  turbulence  statistics 
over  an  area  that  is  much  greater  than  the  spacings  between  the  individual  roughness  ele¬ 
ments  comprising  the  urban  canopy.  This  provides  a  simplified  one-dimensional  mean  wind 
and  turbulence  model  for  urban  canopy  flow  whereby  the  flow  statistics  depend  only  on  the 
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height  2.  In  this  simplified  description  of  urban  flow,  the  spatially-averaged  mean  wind  and 
turbulence  vary  only  with  height,  and  the  detailed  three-dimensional  description  of  flow 
and  turbulence  characteristics  about  the  individual  buildings  is  lost. 

The  primary  objective  of  this  report  is  to  describe  and  test  a  simple,  analytical  one¬ 
dimensional  model  for  the  prediction  of  mean  flow  and  turbulence  within  and  above  an 
urban  canopy.  This  model  was  developed  originally  for  vegetative  canopies  (e.g.,  forests, 
cereal  crops,  etc.),  but  will  be  applied  here  to  urban  canopies.  While  there  are  similarities 
between  vegetative  and  urban  canopies,  there  are  enough  important  differences  between 
these  two  types  of  canopies  (e.g.,  characteristic  density  of  roughness  elements  and  their 
stiffness  or  ability  to  react  to  the  drag  forces)  to  require  an  independent  evaluation  of 
the  model  for  the  urban  canopy  application.  It  is  hoped  that  this  simple,  analytical  one¬ 
dimensional  model  for  urban  canopy  flows  can  be  used  to  provide  useful  input  information 
on  the  statistics  of  the  wind  field  required  to  drive  LS  models  of  urban  dispersion. 
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Model  for  mean  velocity  field  and  shear  stress 


If  the  terrain  is  reasonably  homogeneous  at  the  mesoscale,  Monin-Obukhov  similarity  theory 
provides  a  simple  and  practical  method  for  the  description  of  the  mean  wind  profile  and 
shear  stress  above  the  roughness  sublayer,  where  the  perturbations  on  the  wind  arising 
from  the  roughness  elements  have  decayed.  In  this  section,  we  are  concerned  with  the 
specification  of  simple,  one-dimensional  models  for  the  spatially-averaged  mean  wind  and 
shear  stress  within  the  urban  canopy  (viz.,  within  the  roughness  elements  comprising  the 
urban  canopy).  Unless  otherwise  indicated,  all  flow  variables  of  interest  in  what  follows 
refer  to  the  bulk  or  spatially-averaged  characteristics  of  the  flow  obtained  by  averaging 
over  a  thin,  horizontal  slab  with  a  characteristic  length  scale  that  is  much  greater  than  the 
characteristic  spacings  between  the  individual  obstacles.  The  urban  canopy  is  assumed  to 
be  horizontally  uniform,  so  that  after  the  spatial  averaging  the  mean  flow  and  turbulence 
within  the  canopy  show  a  systematic  variation  in  the  vertical  direction  only. 

To  begin,  a  particular  ansatz  is  chosen  to  model  the  mean  wind  speed  U(z )  within  the 
canopy;  namely, 

U(z)  =  U(hc)  exp(— a(l  -  z/hc)),  (1) 

where  a  is  the  attenuation  coefficient  and  hc  is  the  height  of  the  canopy.  This  is  the 
exponential  profile  that  was  first  suggested  by  Inoue  [1]  and  Cionco  [2]  as  an  appropriate 
wind  profile  within  vegetative  canopies.  Cionco  [3]  has  obtained  values  of  a  for  different 
vegetative  canopies  by  fitting  the  ansatz  of  Equation  (1)  to  measured  velocity  profiles.  He 
indicated  that  a  lies  in  the  range  between  1  and  3  for  moderately  dense  arrays  of  semi¬ 
rigid  roughness  elements  (such  as  mature  corn  and  some  types  of  trees).  It  should  be 
noted  that  it  is  possible  to  “derive”  the  empirical  exponential  wind  profile  by  starting  with 
the  mean  streamwise  momentum  balance  equation  [cf.  Equation  (2)  below],  using  a  flux 
gradient  model  for  the  shear  stress  [u'w'  =  —KmdU/dz,  where  Km  is  a  turbulent  transport 
coefficient  (eddy  viscosity)],  and  using  the  Prandtl-von  Karman  mixing  length  theory  to 
specify  the  eddy  viscosity  as  Km  =  l2dU/dz,  where  the  mixing  length  l  is  most  simply 
taken  as  a  constant  within  the  canopy. 

Assuming  that  the  flow  is  stationary  with  no  mean  streamwise  pressure  gradient,  the  mean 
momentum  equation  for  the  spatially-averaged  mean  wind  U  within  a  canopy  reduces  to  a 
balance  between  the  vertical  divergence  of  the  kinematic  shear  stress  r(z)  =  —u'w'(z)  and 
the  form  or  pressure  drag,  so 

J =CdAU2(z ),  (2) 
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where  Cd  is  the  obstacle  in  situ  drag  coefficient  that  accounts  for  the  effects  of  the  drag  of 
canopy  turbulence  and  also  for  the  sheltering  within  the  canopy,  and  A  is  the  frontal  area 
density  which  is  defined  as  the  frontal  area  of  the  individual  obstacles  per  unit  volume  of 
the  canopy.  In  view  of  Equations  (1)  and  (2),  it  follows  that  the  shear  stress  must  have  the 
following  form  since  its  ^-derivative  is  proportional  to  the  square  of  the  mean  wind  speed: 

-u'w'(z)  =  u*exp(— 2a(l  -  z/hc)),  (3) 

where  u*  =  (~u'w'(z  >  hc))1,/2  is  the  friction  velocity  above  the  canopy.  It  is  assumed 
implicitly  that  a  constant  stress  layer  exists  above  the  canopy.  If  we  substitute  Equation  (3) 
into  Equation  (2) ,  it  is  found  that  the  attenuation  coefficient  has  the  explicit  form 

_  CdAhc 

~  2ui/u*{hcy  (4) 

It  is  worth  noting  that  a  is  the  ratio  of  the  non-dimensional  drag  parameter  CdAhc  (which 
embodies  information  about  the  bulk  canopy  structure)  and  the  skin  friction  2u2/U2(hc) 
(which  contains  information  about  the  interaction  of  the  airflow  with  the  canopy  at  canopy 
height  and,  as  such,  can  be  interpreted  as  a  bulk  drag  coefficient). 

Massman  [4]  showed  how  the  model  for  the  mean  wind  U  and  shear  stress  r  can  be  gener¬ 
alized  to  accomodate  an  arbitrary  canopy  structure,  viz.,  a  canopy  where  Cd  and/or  A  can 
be  functions  of  z) .  The  modification  of  the  mean  wind  and  shear  stress  profiles  in  this  case 
is  minor — it  merely  involves  substituting  for  z/hc  in  Equations  (1)  and  (3),  the  variable 
((z)/((hc).  Here,  ((z)  is  taken  as  the  new  independent  variable  that  is  defined  in  terms  of 
height  2  above  the  ground  surface  beneath  the  canopy  elements  as 

C(z)  =  T  Cd(z')A(z')dz'.  (5) 

Jo 

Furthermore,  in  the  case  of  a  canopy  of  arbitrary  structure,  the  attenuation  coefficient  in 
Equation  (4)  has  ((hc)  substituted  for  the  non-dimensional  drag  parameter  CdAhc  in  the 
numerator.  However,  because  such  detailed  information  on  Cd(z)  and/or  A{z)  is  usually 
unavailable,  we  prefer  to  represent  the  bulk  canopy  structure  using  a  constant  effective  or 
bulk  drag  parameter  ( CdA)b ,  which  for  the  case  of  a  canopy  with  arbitrary  structure  Cd(z) 
and  A(z)  would  simply  be  defined  as  [cf.  Equation  (2)] 

T(hc»u!  =  i  =  (C^y  (£)’  d±  ,  /‘ (£)*  -i.  (6) 

In  fact,  in  all  cases  Cd (z)A(z)  is  simply  treated  as  an  adjustable  model  function  (rather 
than  as  a  strict  model  input)  that  is  obtained  by  tuning  the  function  to  match  either  the 
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model  prediction  of  the  mean  wind  against  available  observations,  or  by  diagnosing  the 
function  from  the  observed  shear  stress  divergence.  Currently,  there  does  not  exist  a  theory 
that  could  be  used  to  predict  the  drag  parameter  from  a  knowledge  of  the  canopy  topology 
and  the  aerodynamic  drag  of  the  roughness  elements  (obstacles)  in  isolation.  In  view  of 
these  uncertainties,  there  does  not  seem  to  be  much  gained  in  trying  to  use  a  local  drag 
parameter  Cd{z)A{z)  in  the  flow  specification  in  place  of  the  simpler  constant  effective  drag 
parameter  CdA. 

The  specification  of  the  mean  wind  profile  of  Equation  (1)  or  the  shear  stress  profile  of 
Equation  (3)  can  be  used  to  determine  the  displacement  height  d  for  the  canopy.  Thom 
[5]  first  observed,  and  Jackson  [6]  subsequently  demonstrated  theoretically  that  d  can  be 
interpreted  as  the  mean  height  of  momentum  absorption  by  the  canopy,  which  in  view  of 
Equation  (2)  can  be  also  interpreted  as  the  level  at  which  the  mean  drag  force  in  the  canopy 
appears  to  act.  Hence, 

d  f1  —u'w'  z 

hc  J0  ul  hc 


Insertion  of  either  Equation  (3)  into  the  first  part  of  Equation  (7)  or,  equivalently,  insertion 
of  Equation  (1)  into  the  second  part  of  Equation  (7)  leads  to  the  following  explicit  expression 
for  the  displacement  length: 

d/hc  =  1  -  ^  (l  -  exp(-2a)).  (8) 

Finally,  the  model  here  can  be  used  to  provide  an  estimate  for  the  roughness  length  of  the 
canopy  in  terms  of  the  attenuation  coefficient  a  and  the  drag  parameter  CdA.  If  we  assume 
that  the  semi-logarithmic  mean  wind  profile  in  the  inertial  sublayer  can  be  extrapolated 
down  to  the  top  of  the  canopy  at  z  =  hc,  we  have 

z0/hc  =  (1  -  d/hc)  exp  >  (9) 

where  k  is  von  Karman’s  constant.  Now,  if  we  substitute  for  d/hc  from  Equation  (8)  and 
I7(hc)/ti*  from  Equation  (4),  we  get 

zo/hc  -  ^(l-exp(-2a))exp(-K(2a/CdAhc)1/2y  (10) 
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Model  for  the  TKE  and  velocity  variances 


The  key  to  the  determination  of  the  velocity  variances  centers  on  the  construction  of  an  an¬ 
alytical  solution  to  a  simplified  transport  equation  for  the  turbulence  kinetic  energy  (TKE) 
k.  Massman  and  Weil  [7]  showed  how  to  construct  a  simple  analytical  model  for  k  based 
on  the  second-order  closure  model  for  canopy  flows  developed  by  Wilson  and  Shaw  [8].  The 
transport  equation  for  k  =  (erf  +  <rf  +  erf,)/ 2  (where  <JU1  av,  and  aw  are  the  streamwise  (x), 
cross-stream  (y),  and  vertical  (z)  components  of  the  velocity  standard  deviations,  respec¬ 
tively)  implied  by  the  second-order  closure  model  for  turbulence  statistics  of  Wilson  and 
Shaw  [8]  assumes  the  following  form: 


d_ 

dz 


w + H  w  -  <*/»)+**"• + -  °’ 


(10) 


where  erf  is  twice  the  turbulence  kinetic  energy  (i.e.,  k  =  cf/2).  From  left  to  right,  the  terms 
in  the  transport  equation  for  erf  represent  the  following  physical  processes:  (1)  turbulent 
transport  of  erf;  (2)  rate  of  viscous  dissipation  of  erf,  which  converts  turbulence  kinetic 
energy  to  heat;  (3)  return-to-isotropy,  expressing  the  rate  at  which  turbulence  energy  is 
redistributed  between  the  three  components  of  the  velocity  variance;  (4)  wake  production, 
accounting  for  the  production  of  turbulence  kinetic  energy  from  the  mean  kinetic  energy 
of  the  flow  in  the  wakes  of  the  roughness  elements  (obstacles);  and,  (5)  shear  production, 
accounting  for  the  conversion  of  the  mean  kinetic  energy  to  turbulence  kinetic  energy  due 
to  the  action  of  the  shear  stress  against  the  mean  velocity  shear.  It  should  be  emphasized 
that  the  return-to-isotropy  term  does  not  normally  appear  in  the  transport  equation  for  the 
turbulence  kinetic  energy,  but  its  presence  in  this  case  results  from  the  particular  closure 
assumptions  adopted  by  Wilson  and  Shaw  [7]  for  their  second-order  closure  model  for  canopy 
flow. 


To  close  Equation  (10),  it  is  necessary  to  parameterize  the  three  canopy  length  scales  A* 
(i  =  1, 2, 3).  Wilson  and  Shaw  [7]  made  the  assumption  that  these  three  length  scales  vary 
with  height  z  within  the  canopy  in  a  similar  fashion,  so  A^z)  =  V{l(z ),  where  Vi  ( i  =  1, 2, 3) 
are  constants  and  l  is  a  characteristic  length  scale.  Generally,  the  parameterization  scheme 
used  by  Wilson  and  Shaw  prescribed  the  characteristic  length  scale  l(z )  as  a  function  of  the 
density  and  drag  properties  of  the  (vegetative)  canopy  and  height  above  the  surface  such  that 
the  length  l  is  continuous  and  assumes  the  maximum  possible  value  allowed  by  the  following 
constraints:  (1)  l  =  0  at  z  =  0;  (2)  l  <  a+/(CdA)  for  0  <  z  <  hc‘,  and,  (3)  \dl/dz\  <  k  for 
z  >  0,  where  a*  is  a  constant,  with  a*  and  Cd  chosen  to  obtain  the  best  fit  to  observed 
canopy  flow  data.  However,  to  ensure  that  an  analytical  solution  can  be  constructed  for 
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Equation  (10),  Massman  and  Weil  [7]  relaxed  some  of  these  constraints  and  simply  assumed 
that  the  characteristic  length  scale  can  be  parameterized  as  l(z)  =  a*/ ( Cd(z)A(z )).  Note 
that  for  the  case  we  are  interested  in,  l(z)  is  a  constant  throughout  the  canopy  since  C'dA  is 
assumed  to  be  a  constant.  The  constants  z/j  in  the  definition  of  Xt  are  chosen  to  be  consistent 
with  the  constant  stress  layer  and  semi-logarithmic  mean  wind  profile  that  is  assumed  to 
exist  immediately  above  the  canopy.  In  particular,  it  is  assumed  that  immediately  above 
the  canopy,  the  velocity  standard  deviations  cq  are  constant  and  proportional  to  the  friction 
velocity  (i.e.,  al  —  jiU *,  where  7 »  (i  =  1,2, 3)  are  constants  that  are  assumed  to  be  known 
from  the  properties  of  the  constant  stress  layer  above  the  canopy).  To  be  consistent  with  the 
specified  properties  of  the  constant  stress  layer,  the  closure  coefficients  17  must  be  chosen 
as  follows:  vx  =  (7 1  +  7!  +  7a)-172.  ^3  =  (Ti  +  I2  +  7 !)3/2>  and  v2  =  ^3/6  -  72/(2i/i)- 
With  these  relationships  for  the  closure  coefficients  17  and  the  parameterization  for  the 
characteristic  length  scale  /,  the  transport  equation  for  ae  in  Equation  (10)  reduces  to 
(assuming  a  constant  drag  parameter  CdA ) 


— 1  =  A \CdAfc\  -  — ( CdAfU 3 


atux  dz 


(11) 


dz2  '  '  ~  a»z/i 

where  A2  =  ‘iu2 j a2.  Massman  and  Weil  [7]  showed  that  Equation  (11)  possesses  the 
following  analytical  solution  subject  to  the  upper  boundary  condition  at  z  =  hc  specified 
by  (rl{z  =  he)/ul  =  v%\ 
k1/2  _  1  ae 
u*  y/2  u * 

_  _1_  l/3e-AC’dAZic(l-z//ic)  _|_  ^  /'e-3a(l-z//lc)  _  e-ACdAhc(l-z/hc)\  173  ^  ^2) 

where  =  —  (9it*/C/(/ic))/(2a*^i[9/4  -  A2tz*/C/4(/ic)]). 


The  expression  for  the  turbulence  kinetic  energy  k  in  Equation  (12)  can  be  used  to  obtain 
predictions  for  the  velocity  variances.  The  simplest  method  is  simply  to  derive  the  velocity 
variances  from  the  calculated  turbulence  kinetic  energy  by  assuming  equilibrium  partition¬ 
ing.  In  this  approach,  it  is  assumed  that  the  component  velocity  variances  in  the  canopy 
are  proportional  to  the  turbulence  kinetic  energy  and  that  the  proportionality  constants 
are  the  same  as  those  in  the  constant  stress  layer  overlying  the  canopy.  In  consequence, 
equilibrium  partitioning  provides  the  following  predictions  for  the  standard  deviations  of 
the  velocity  within  the  canopy: 

01(2)  =  J~2  —  ^2  +  j2y/2  (2^(z))172’  z  -  he,  (13) 

where  k  is  determined  by  Equation  (12).  Note  that  using  equilibrium  partitioning,  the 
turbulent  velocity  variances  are  scaled  versions  of  the  turbulence  kinetic  energy  and,  hence, 
have  exactly  the  same  shape. 
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An  alternative  to  equilibrium  partitioning  for  the  derivation  of  the  velocity  variances  would 
be  to  apply  an  algebraic  stress  model  (ASM).  The  algebraic  stress  model  that  is  used  here 
is  derived  from  the  well  understood  and  successful  second-order  closure  model  developed 
by  Launder,  Reece,  and  Rodi  [9]  with  the  incorporation  of  additional  terms  to  describe  the 
physics  particular  to  canopy  flows.  For  three-dimensional  canopy  flows,  it  can  be  shown 
after  some  rather  lengthy  and  tedious  algebraic  manipulations  (Yee,  1999,  unpublished) 
that  this  ASM  leads  to  the  following  form  for  the  anisotropy  tensor  a,ij  =  u^u'Jk  —  |<5» j, 
(Sij  is  the  Kronecker  delta  function): 

“  "b  ( Sik&kj  T  SjkQ/ki  —SlkClkl5ij  \  “I”  63  ( T  jk&ki)  ^1  (14) 


where 


=  I*  (dUi 

~  2  e  \  dxj 


is  the  normalized  mean  rate-of-strain  tensor  (viz.,  normalized  by  the  turbulence  time  scale 
fc/e,  where  e  is  the  viscous  dissipation  rate); 


13  2  e  \  &Xj  dxi  ) 


is  the  normalized  mean  vorticity  (rotation)  tensor; 


Ci  +  Ps/ e  +  Pw/ e  —  1 


.  8  t  5-9. 

=  is’  ^““ir 


7c2+  1 


are  numerical  constants.  Furthermore,  ci  and  C2  are  closure  constants  that  arise  in  the 
Launder,  Reece,  and  Rodi  [9]  parameterization  of  the  pressure-strain  term.  Note  that  g* 
[cf.  Equation  (17)]  embodies  the  information  on  the  net  production  (shear  Ps  and  wake 
Pw )  per  unit  dissipation  of  the  TKE.  Equation  (14)  constitutes  a  closed  system  of  algebraic 
equations  for  the  determination  of  the  Reynolds  stress  anisotropy  in  terms  of  the  mean  ve¬ 
locity  gradients.  Once  the  mean  flow  and  TKE  and  its  dissipation  are  known,  Equation  (14) 
can  then  be  solved  numerically  at  each  point  in  the  flow  domain  for  the  anisotropy  stress 
components  (whose  normal  stress  components  along  the  diagonal,  in  turn,  determines  the 
partitioning  of  the  TKE). 


Equation  (14)  simplifies  considerably  for  the  spatially-averaged  flow  variables  (mean  wind 
and  turbulence  statistics)  describing  an  effectively  horizontally  homogeneous  canopy.  For 
flow  variables  that  depend  only  on  z,  Equations  (14)  to  (18)  reduce  to  the  following  form: 


8  +  12c2 


x  (— 25oOi3); 
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(20) 


<222  I8C2  —  10  /  „  x 

—  =  ^3—  X  (-2S„a13); 

af  =  x  (-2 Scan).  (21) 

Here,  an  =  <72/fc-2/3,  022  =  cr^/ fc— 2/3,  and  <233  =  cr2/fc-2/3.  Also,  So  =  \\dU  fdz  is  the 
normalized  mean  rate-of-strain,  ai3  =  u'w'/k  is  the  normalized  kinematic  shear  stress,  and 
g*  is  given  by  Equation  (17)  with  Ps  =  —u'w’dU/dz  (shear  production)  and  Pw  =  ^ CdAU 3 
(wake  production,  or  work  done  by  the  mean  flow  against  form  drag).  To  complete  the  ASM 
partitioning  of  the  velocity  variances,  we  need  a  model  for  the  TKE  dissipation  rate  e(z). 
Although  Massman  and  Weil  [7]  provide  a  model  for  e(z),  we  prefer  to  use  a  TKE  dissipation 
rate  model  suggested  by  Wilson,  Finnigan,  and  Raupach  [10]  that  explicitly  distinguishes 
between  shear-generated  and  wake-generated  turbulence.  Here,  the  TKE  dissipation,  e,  is 
modelled  as  e  =  max(ecc,  efd)  where: 

(c  Ad3/2 

ecc  -  ~-A;  €fd  =  P*CdAUk .  (22) 

The  term  ecc  is  the  standard  parameterization  for  free-air  viscous  dissipation,  the  only  TKE 
sink  above  the  canopy.  The  specification  for  efd  (‘form  drag’),  which  converts  resolved  TKE 
to  the  small,  rapidly  dissipated  ‘wake-scales’,  can  be  justified  on  physical  grounds.  The 
constant  ce  =  ul/k(z  —  hc)  is  the  shear  stress  to  TKE  ratio  at  the  canopy  top  (z  =  hc). 
The  closure  coefficient  /?*  =  1  was  suggested  by  Wilson,  Finnigan,  and  Raupach  [10].  The 
turbulence  length  scale  A  is  determined  from  A  =  max(Af,  A0),  where 


1_  _J_  J_  1  1  1 

A i  kz  Ac  ’  A0  k(z  —  d)  Loo  ’ 


(23) 


where 

Ac  =  (k(z  -  ^c))1/2/ dU/dz\z=hc  (24) 

is  a  canopy  shear  length  scale  measuring  the  strength  of  the  wind  shear  at  the  canopy  top 
that  Raupach,  Finnigan,  and  Brunet  [11]  have  identified  as  being  a  critical  length  scale 
governing  the  eddy  transport  within  a  canopy.  Here,  \  is  an  ‘inner’  scale,  Ac  is  an  ‘outer’ 
scale,  and  L0 Q  is  an  imposed  limiting  length  scale  that  restricts  the  depth  of  the  inertial 
sublayer  overlying  the  canopy. 


The  closure  coefficients  Ci  and  C2  arising,  respectively,  from  the  return-to-isotropy  and 
rapid  part  of  the  Launder,  Reece,  and  Rodi  [9]  parameterization  of  the  pressure-strain 
redistribution  term  in  the  transport  equation  for  the  Reynolds  stress  tensor,  can  be  set  by 
the  physical  properties  of  the  flow  in  the  logarithmic  constant  stress  layer  assumed  to  overlie 
the  canopy.  Again,  if  in  this  constant  stress  layer,  erf  =  (i  —  1,2,3  =  u,v,w),  then  it 
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follows  that  (assuming  in  this  layer  that  the  production  of  TKE  is  equal  to  its  dissipation, 
corresponding  to  the  condition  of  local  equilibrium) 


Here,  c~ 1  =  ^(7^  +  7^  +  7^)  is  the  prescribed  (observed)  equilibrium  value  at  canopy 
top  of  the  ratio  of  the  TKE  to  the  shear  stress  [cf.  Equation  (22)].  Given  the  relations  in 
Equations  (25)  to  (27),  specified  values  for  7 u,  jv,  and  7™  will  yield  values  for  ci,  C2,  and  ce. 
Typically,  for  most  turbulent  boundary  layers,  ju  ~  2.4,  jv  ~  1.9,  and  jw  ~  1.25  implying 
that  ce  ta  0.18.  For  a  wide  range  of  thin  shear  flows  in  equilibrium,  the  shear  stress/TKE 
ratio  (i.e.,  the  constant  ce )  is  relatively  constant.  Indeed,  various  measurements  which  have 
been  summarized  by  Townsend  [12]  indicate  that  for  boundary  layers,  wakes,  and  mixing 
layers,  the  ratio  is  nearly  the  same  and  varies  between  about  0.2  to  0.3.  The  constant  ce 
in  the  engineering  literature  is  often  referred  to  as  Bradshaw’s  constant,  and  sometimes  as 
Townsend’s  constant  (although  the  latter  is  defined  in  our  current  notation  to  be  cej 2). 
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Parameterizations  for  model  inputs 


The  key  inputs  to  the  model  are  the  “skin  friction”,  u2/U2(hc),  (actually,  the  skin  friction 
Cf  is  twice  this  value)  and  either  the  attenuation  coefficient,  a,  or  the  normalized  drag 
parameter,  CdAhc  [cf.  Equation  (4)]  Furthermore,  we  need  to  know  the  equilibrium  values 
of  k/u 2  (which  we  have  defined  to  be  cj1),  and  cq/u*  (i  —  1, 2, 3  =  u,  v,  w)  at  the  canopy 
top.  Only  two  of  the  values  cq/u*  need  to  be  known  at  the  canopy  top  in  addition  to  k/u2 
since  k/u2  =  \{cr2u/u\  +  cr2/ul  +  (t2/u2).  However,  the  critical  input  values  are  u2jU2(hc) 
and  CdAhc  (or,  a),  since  the  values  of  cq/u*  (i  =  1, 2, 3)  and/or  ce  does  not  appear  to  vary 
greatly  from  canopy  to  canopy,  or  from  one  local  equilibrium  layer  to  another. 


If  observed  values  for  ul/U2(hc )  and  CdAhc  (or,  a)  are  available  for  a  particular  canopy, 
then  they  can  be  used  directly  in  the  model.  For  general  applications  of  the  model,  when 
observed  values  for  u2/U2(hc)  and  a  are  not  available,  parameterizations  for  these  quantities 
are  required.  For  vegetative  canopies,  Massman  [4]  suggested  the  following  parameterization 
for  the  skin  friction  at  canopy  top: 


u*/U(hc)  =  ki  —  k2exp(—k3CdAhc), 


(28) 


where  k\  =  0.320,  k2  —  0.264,  and  k3  =  15.1  are  model  constants.  Also,  Raupauch  [13,  14] 
developed  for  random  arrangements  of  obstacles  that  are  not  too  densely  spaced  together,  a 
parameterization  for  the  skin  friction  that  is  expressed  as  an  implicit  function  of  the  frontal 
area  index  \F  (i.e.,  ratio  of  the  frontal  area  of  the  elements  facing  the  wind  to  the  total 
ground  area  occupied  by  the  elements): 


U(hc)/u* 


exp(c\FU(hc)  /  (2u*)) 

(cs  +  ±cd\F)1/2 


(29) 


where  Cs  is  the  drag  coefficient  of  the  substrate,  Cd  is  the  drag  coefficient  of  the  obstacles 
in  isolation,  and  c  «  0.5  is  an  empirical  shelter  coefficient  that  is  determined  by  the  rate  at 
which  a  roughness  element  wake  spreads  in  the  cross-stream  direction.  Because  Raupach’s 
model  cannot  properly  describe  overlapping  of  sheltering  in  densely  packed  roughness  arrays, 
it  can  only  be  used  for  XF  less  than  about  0.1  to  0.2. 


For  regularly  spaced  arrays  of  bluff  roughness  elements  (urban  canopies) ,  much  of  the  effort 
has  focussed  on  the  parameterization  of  the  roughness  length,  zo,  as  a  function  of  the 
frontal  area  index,  \F.  However,  these  parameterizations  for  Zo  must  imply  implicitly 
also  a  parameterization  for  the  skin  friction.  In  particular,  all  these  parameterizations  start 
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with  the  assumption  that  the  semi-logarithmic  velocity  profile  in  the  inertial  sublayer  (lying 
above  the  canopy)  can  be  extrapolated  down  to  the  canopy  top  at  z  =  hc  such  that  the 
skin  friction  u2/U2(hc)  is  linked  to  z0  as 


£o 

hc 


(30) 


The  error  in  this  assumption  is  small,  and  is  discussed  in  detail  in  Hanna  and  Britter  [15]. 
Many  parameterizations  for  zq  involve  modeling  U{hc)/u *  as  a  function  of  the  topological 
parameters  that  characterize  the  urban  canopy.  For  example,  the  Macdonald  et  al.  [16] 
parameterization  for  z0  for  regular  arrays  of  obstacles  (square  and  staggered)  effectively 
uses  the  following  expression  for  the  skin  friction  in  Equation  (30): 

= \^Cd  (/  -  £)  Af’  <si> 


where  Cd  is  the  drag  coefficient  of  an  isolated  roughness  element  mounted  on  the  surface 
and  (3  is  a  correction  factor  (empirical  parameter)  for  the  drag  coefficient  that  includes  the 
influences  of  the  incident  turbulence  length  scale,  the  incident  wind  angle,  the  presence 
of  rounded  corners  in  the  obstacles,  and  the  geometry  of  the  roughness  array.  For  square 
and  staggered  arrays  of  cubical  obstacles,  Macdonald  et  al.  [17]  calibrated  this  relationship 
using  their  wind  tunnel  data  and  recommended  values  of  0.55  and  1.0  for  (3,  for  square 
and  staggered  arrays  of  cubical  obstacles,  respectively.  An  alternative  parameterization  has 
been  suggested  by  Duijm  [16]  who  uses 


<  1 


U2{hc)  2CdXF’eff' 


(32) 


in  the  parameterization  for  z0/hc  in  Equation  (30).  Here,  A  p,ef{  is  an  effective  frontal  area 
index  that  accounts  for  the  sheltering  and  assumes  the  form 


^F,eff  —  Af  ^1  -  >  (33) 

where  S  is  an  empirical  sheltering  efficiency  factor  of  order  one  and  dpi  is  an  in-plane 
displacement  height  that  was  originally  introduced  by  Bottema  [18,  19]  to  account  for  the 
effect  of  sheltering  and  the  presence  of  recirculation  zones. 


It  is  not  possible  based  on  the  available  data  to  recommend  a  particular  parameterization 
for  zq.  The  large  uncertainties  in  the  determination  of  zo  (especially  from  mean  velocity 
profiles  above  the  canopy)  and  the  large  scatter  in  the  measurements  of  zo  (c.f.  Grimmond 
and  Oke  [20])  do  not  provide  a  standard  that  can  be  used  to  distinguish  between  the 
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different  parameterizations  for  Zo  proposed  in  the  literature.  Indeed,  Bottema  [18]  and 
Grimmond  and  Oke  [20]  provide  a  comprehensive  survey  of  data  on  zo  obtained  from  various 
urban  canopies,  and  strongly  emphasize  that  there  are  numerous  reasons  for  doubting  their 
correctness.  In  consequence,  the  various  parameterizations  for  zo  represent  nothing  more 
than  ad-hoc  formulae  that  can  only  really  be  justified  by  their  success  in  fitting  the  limited 
observations  used  for  their  development. 

The  remaining  parameter  that  needs  to  be  specified  is  the  attenuation  coefficient  a.  Best 
values  of  a  for  numerous  canopies  (but,  primarily  vegetative)  have  been  compiled  by  Cionco 
[3].  More  recently,  Macdonald  et  al.  [21]  found  least-squares  fitted  values  for  a  for  regular 
roughness  arrays  (urban  canopies)  of  cubical  obstacles  with  plan  area  indices  (i.e.,  ratio  of 
obstacle  plan  area  to  total  ground  area  occupied  by  the  obstacles)  in  the  range  from  about 
0.05  to  0.30.  These  researchers  used  this  data  to  recommend  the  following  relationship 
between  a  and  the  plan  area  index  Ap  for  Ap  £  (0.05, 0.30): 

a  =  -0.24  +  12.5Ap.  (34) 
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Results 


Massman  and  Weil  [7]  tested  the  one-dimensional  analytical  model  for  canopy  flow  against 
a  number  of  horizontally  homogeneous  vegetative  canopies,  primarily  to  determine  an  ap¬ 
propriate  value  for  the  closure  coefficient  a*  required  to  define  the  characteristic  length 
scale  l  (z) .  To  this  purpose,  these  researchers  calibrated  the  model  against  corn  canopy  data 
from  Shaw  et  al.  [22]  and  Wilson  et  al.  [23],  spruce  and  aspen  forest  data  from  Amiro 
[24],  Douglas  fir  forest  canopy  data  from  Lee  and  Black  [25],  pine  forest  data  from  Katul 
and  Albertson  [26],  and  a  wind  tunnel  ‘wheat’  canopy  data  from  Brunet  et  al.  [27].  This 
calibration  procedure  involved  some  rather  fine-tuning  to  the  data  consisting  of  adjusting 
the  variable  drag  coefficient  Cd(z )  to  provide  the  best  fits  of  the  model  predictions  to  the 
observed  canopy  profiles  of  mean  wind  U(z)  and  shear  stress  —u'w'(z),  and  then  adjusting 
the  parameter  a*  to  provide  an  optimal  fit  to  the  observed  canopy  profiles  of  the  turbulent 
velocity  standard  deviations  au  and  aw.  Based  on  this  trial  and  error  adjustment  to  the 
observed  data,  Massman  and  Weil  [7]  found  that  the  values  for  a*  were  generally  between 
about  0.04  and  0.06.  Using  these  results,  they  recommended  that  a  value  of  a«  =  0.05  be 
used  for  general  model  application. 

In  this  section,  we  make  model  comparisons  with  data  from  “urban”  canopies  (i.e.,  roughness 
element  arrays  composed  of  solid  bluff  obstacles,  rather  than  flexible  elements  that  are 
characteristic  of  vegetative  canopies).  Unfortunately,  there  is  a  dearth  of  data  on  observed 
profiles  of  mean  wind  speed  and  turbulence  statistics  within  an  urban  canopy.  The  first 
urban  canopy  considered  is  described  in  detail  by  Raupach,  Coppin,  and  Legg  [28].  This 
canopy  (referred  to  as  the  Tombstone  canopy)  consisted  of  an  extensive  array  of  thin, 
“billboard-shaped”  obstacles  (or,  vertical  bars)  with  dimensions  60  mm  high,  10  mm  wide, 
and  1  mm  thick,  arranged  in  a  staggered  pattern  with  a  spacing  of  60  mm  in  the  cross- 
stream  direction  and  44  mm  in  the  streamwise  direction  to  give  a  frontal  area  index  for  the 
array  of  0.23.  For  specific  comparisons  between  the  model  predictions  and  the  observations, 
the  drag  parameter  CdAhc,  skin  friction  ut/U(hc).  and  the  7,’s  (i.e.,  ratio  of  the  velocity 
standard  deviation  in  the  i-th  direction  to  the  friction  velocity  at  the  canopy  top)  were 
taken  from  the  observations. 

Figures  1  and  2  compare  measured  profiles  of  the  mean  wind  speed  and  kinematic  shear 
stress  for  the  Tombstone  canopy  with  the  model  predictions.  It  is  seen  that  the  attenuation 
coefficient,  a,  inferred  from  the  observed  values  of  the  drag  parameter  and  skin  friction 
factor  [cf.  Equation  (4)]  is  too  large,  causing  the  predicted  mean  wind  speed  and  shear 
stress  to  decay  too  rapidly  within  the  canopy.  Least-squares  fitting  of  the  exponential  profile 
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(cf.  Equation  (1))  to  the  observed  profile  of  U(z )  gave  a  value  of  1.41  for  a,  rather  than  the 
larger  inferred  value  of  1.78.  Using  this  fitted  value  for  a  gave  reasonable  predictions  for 
the  mean  wind  speed  and  shear  stress  within  the  Tombstone  canopy.  In  consequence,  for 
the  comparisons  of  model  predictions  with  observations  for  the  turbulence  kinetic  energy 
and  turbulent  velocity  standard  deviations,  we  will  use  a  =  1.41  in  the  analytical  model 
unless  otherwise  indicated. 

Figure  3  shows  a  comparison  of  the  observed  and  modeled  turbulence  kinetic  energy  for  the 
Tombstone  canopy.  Using  the  value  of  a*  =  0.05  recommended  by  Massman  and  Weil  [7],  it 
is  seen  that  the  model  predictions  for  k  with  attenuation  coefficients  of  1.78  (inferred  value) 
and  1.41  (least-squares  fitted  value)  provide,  respectively,  a  lower  and  upper  bound  for  the 
observed  values  of  k.  However,  the  model  profile  for  k  with  a*  =  0.055  (chosen  by  trial 
and  error)  and  the  least-squares  fitted  value  of  a  =  1.41  seems  to  provide  the  best  visual 
comparison  with  the  observed  k  profile.  However,  even  so,  this  model  profile  underestimates 
k  in  the  region  at  about  mid-canopy  height. 

Figures  4  and  5  show  the  measured  and  modeled  profiles  for  au  and  am,  respectively.  The 
modeled  profiles  for  the  normal  stresses  here  were  obtained  using  equilibrium  partitioning 
of  the  turbulence  kinetic  energy.  It  is  seen  that  agreement  between  observation  and  model 
predictions  for  au  for  the  case  with  a  =  1.41  and  a*  =  0.055  is  good.  However,  the  model 
profile  for  aw  for  this  case  does  not  agree  as  well  with  the  observed  profile,  primarily  because 
the  observed  profile  for  aw  shows  more  concavity  than  that  for  au.  In  particular,  the  model 
profile  for  aw  with  a  =  1.41  and  a*  =  0.055  underestimates  the  observed  values  in  the  upper 
part  of  the  canopy.  Because  the  model  profiles  for  the  turbulent  velocity  standard  deviations 
obtained  by  equilibrium  partitioning  are  essentially  constructed  by  applying  scaling  factors 
to  the  associated  model  profile  for  k ,  all  these  profiles  have  necessarily  the  same  shape. 
Hence,  because  of  the  differences  in  shape  of  the  observed  profiles  for  au  and  aw,  it  is  not 
possible  for  the  equilibrium  partitioning  of  k  to  provide  model  profiles  for  the  turbulent 
velocity  standard  deviations  that  can  fit  both  the  observed  profiles  for  au  and  crw. 

Figures  6  and  7  display  the  measured  and  modeled  profiles  for  au  and  <rTO,  respectively.  The 
modeled  profiles  here  are  obtained  by  the  algebraic  stress  partitioning  of  the  TKE.  The 
agreement  of  the  model  and  observations  is  excellent  for  cru  with  a  =  1.41  and  a*  =  0.055, 
but  the  corresponding  model  profile  for  aw  underestimates  the  measurements  for  the  vertical 
velocity  standard  deviation  in  the  middle  to  upper  part  of  the  canopy.  In  spite  of  the  fact 
that  algebraic  stress  partitioning  of  the  TKE  can  produce  different  shapes  for  the  component 
velocity  variance  profiles,  the  agreement  obtained  for  aw  was  no  better  than  that  procured 
from  the  equilibrium  partitioning  (at  least  for  the  Tombstone  canopy). 
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Macdonald  et  al.  [21]  made  measurements  of  the  mean  wind  speed  and  streamwise  turbu¬ 
lent  velocity  standard  deviation  in  urban  canopies  comprised  of  arrays  of  regularly  spaced 
cubical  obstacles  with  plan  area  indices  ranging  from  about  0.05  to  0.9.  These  researchers 
demonstrated  that  for  plan  area  indices  less  than  about  0.2,  the  measured  mean  wind  speed 
within  the  canopy  is  well  modeled  using  the  exponential  wind  profile  of  Equation  (1).  Al¬ 
though  the  plan  area  index  Ap  is  used  to  denote  a  proportion,  the  most  convenient  symbolism 
for  theoretical  interpretation,  most  users  of  this  quantity  prefer  to  do  their  arithmetic  in 
terms  of  percentages,  and  in  calculations  the  same  symbols  will  be  used  even  though  100AP 
would  be  more  correct.  No  confusion  should  be  caused  by  this  slight  ambiguity  of  symbol¬ 
ism.  Figure  8  displays  the  variation  of  the  attenuation  coefficient  a  as  a  function  of  the 
plan  area  index  Ap  for  both  square  and  staggered  arrays,  indicating  that  there  are  only 
minor  differences  between  the  two  array  geometries.  Macdonald  et  al.  [21]  recommended 
the  following  relationship  between  a  and  Ap:  a  =  —0.24  +  0.125AP  (with  Ap  expressed  as  a 
percentage).  However,  this  relationship  predicts  that  a  <  0  for  Ap  <  2  percent.  To  remove 
this  physical  inconsistency,  we  re-fitted  the  data  using  a  power-law  relationship  of  the  form 
a  —  A\p  where  s  is  the  power-law  exponent.  The  best  power-law  fit  to  the  data  is  shown  in 
Figure  8  and  yielded  the  following  result:  a  =  0.104AP  (with  Ap  expressed  as  a  percentage). 
Note  that  in  this  power-law  fit,  s  =  1  almost  exactly. 

The  observed  values  of  the  skin  friction  factor  at  canopy  top  for  the  staggered  and  square 
arrays  are  shown  in  Figure  9  as  a  function  of  the  plan  area  index.  For  each  array  geometry, 
the  skin  friction  is  seen  to  increase  monotonically  for  low  values  of  the  obstacle  density, 
corresponding  roughly  to  the  isolated  flow  regime.  The  skin  friction  increases  with  Ap  to 
a  broad  maximum  that  appears  to  occur  for  Ap  between  about  20  and  30  percent.  This 
range  for  Ap  corresponds  roughly  to  the  wake  interference  regime.  Beyond  Ap  ~  30  percent, 
the  skin  friction  decreases  and  appears  to  approach  an  approximately  constant  value  for 
Ap  >  40  percent,  a  range  which  can  be  associated  with  the  skimming  flow  regime  whereby 
the  obstacles  are  so  closely  packed  that  flow  aloft  does  not  penetrate  into  the  obstacle  array 
and  a  virtually  stagnant  mean  flow  is  present  within  the  array. 

Using  the  values  for  a  and  u*/U(hc )  provided  by  Figures  8  and  9  as  input  to  the  analytical 
model,  Figure  10  shows  the  model  profiles  for  k  for  a  square  array  of  cubicial  obstacles 
with  the  plan  area  index  between  about  0.05  and  0.20.  As  expected,  the  model  predicts  the 
increasing  attenuation  of  k  with  canopy  depth  as  the  obstacle  density  (or,  equivalently,  as 
the  plan  area  index)  increases.  Figures  11  and  12  exhibit  model  predictions  for  au  and  aw, 
respectively,  obtained  by  the  equilibrium  partitioning  of  the  TKE  profiles  in  Figure  10.  As  a 
comparison,  Figures  13  and  14  show  the  model  profiles  for  au  and  aw,  respectively,  obtained 
by  the  algebraic  stress  model  partitioning  of  the  TKE  profiles.The  data  points  in  Figures  11 
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and  13  for  cru/u*  were  extracted  from  Figure  19  of  Macdonald  et  al.  [21].  However,  the 
scale  used  by  Macdonald  et  al.  [21]  for  their  Figure  19  only  permitted  certain  values  of 
cru/u*  at  Ap  =  4.9%  and  20%  to  be  extracted  unambiguously;  in  consequence,  we  have 
shown  only  the  range  in  the  measured  au/u *  data  for  Ap  G  [4.9, 20]%,  as  the  individual  data 
points  at  the  various  plan  area  indices  cannot  be  distinguished  unambiguously  in  Figure  19 
of  Macdonald  et  al.  [21].  Although  the  model  profiles  for  the  turbulent  velocity  standard 
deviations  appear  to  be  similar  for  the  equilibrium  and  algebraic  stress  model  partitioning, 
a  comparison  of  Figures  11  and  13  seems  to  suggest  that  the  latter  partitioning  provides  a 
slightly  better  reproduction  for  the  range  of  variation  in  the  au  observations. 

Macdonald  et  al.  [21]  noted  that  the  exponential  wind  profile  does  not  provide  good  fits 
to  the  observations  of  within  canopy  mean  wind  speed  for  Ap  >  20  percent.  Indeed,  for 
Ap  >  40  percent,  the  obstacle  packing  densities  are  high  enough  to  produce  skimming 
flow  implying  a  near  stagnant  flow  within  the  obstacle  array  with  a  mean  wind  speed 
of  near  zero.  We  do  not  expect  the  analytical  model  to  be  valid  at  these  high  obstacle 
packing  densities.  Nevertheless,  we  applied  the  model  to  obstacle  arrays  with  high  packing 
densities.  A  comparison  of  observed  and  modeled  streamwise  turbulent  velocity  standard 
deviations  for  Ap  =  69%  and  91%  are  displayed  in  Figures  15  and  16,  for  equilibrium  and 
algebraic  stress  model  partitioning  of  the  TKE,  respectively.  Although  the  model  profile 
underpredicts  au  for  Ap  =  69%,  the  prediction  for  Ap  =  91%  surprisingly  seems  to  provide 
a  fairly  good  representation  for  the  observed  profile,  especially  deep  in  the  canopy  where 
(j u  appears  to  be  constant  (approximately  or  better). 
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Conclusions 


In  this  paper,  we  describe  a  fully  analytical  model  for  the  prediction  of  the  one-dimensional 
mean  wind  speed,  kinematic  shear  stress,  turbulence  kinetic  energy,  and  velocity  variances 
in  a  horizontally  homogeneous  canopy.  The  basis  of  the  model  centers  around  an  analytical 
solution  obtained  by  Massman  and  Weil  [7]  for  the  turbulence  kinetic  energy  derived  from 
a  particular  form  of  a  second-order  closure  model  for  canopy  flows  developed  by  Wilson 
and  Shaw  [8].  Using  this  analytical  model  for  the  TKE,  we  showed  how  to  derive  analytical 
expressions  for  the  within  canopy  velocity  variances  using  two  methods:  namely,  (1)  equi¬ 
librium  partitioning;  and,  (2)  algebraic  stress  model  partitioning  derived  from  an  algebraic 
stress  model  version  of  an  existing  second-order  closure  model  for  flows  proposed  by  Laun¬ 
der,  Reece,  and  Rodi  [9]  that  has  been  modified  to  include  the  extra  physics  arising  from 
the  interaction  of  the  airflow  with  the  obstacles  in  the  array. 

The  analytical  model  predictions  have  been  compared  with  some  measurements  of  mean  flow 
and  turbulence  statistics  obtained  in  two  urban  canopies:  namely,  (1)  Tombstone  canopy; 
and,  (2)  arrays  of  regularly  spaced  cubical  obstacles  at  various  plan  area  indices.  It  is  shown 
that  the  model  is  moderately  successful  at  providing  estimates  of  a  number  of  important 
flow  quantities;  namely,  mean  wind  speed  and  various  turbulence  statistics.  Given  the 
limited  data  that  was  available,  we  have  been  unable  to  ascertain  whether  the  algebraic 
stress  model  partitioning  of  the  TKE  offers  any  significant  advantages  in  the  determination 
of  the  velocity  variances  over  the  simpler  equilibrium  partitioning.  Certainly,  the  measured 
profiles  for  au  appear  to  be  different  in  shape  than  those  of  crw  (viz.,  there  appear  to  be 
differences  in  the  concavity  of  the  uu  and  aw  profiles),  a  general  behavior  that  can  be 
captured  potentially  by  the  algebraic  stress  model  partitioning  but  not  by  the  equilibrium 
partitioning  of  the  TKE.  The  assumption  of  an  exponential  wind  profile  may  not  be  valid  for 
urban  arrays  at  high  plan  area  indices  where  the  skimming  flow  regime  is  dominant  and  a 
virtually  stagnant  flow  is  found  within  the  urban  array.  However,  even  under  these  extreme 
conditions,  the  analytical  model  has  been  found  to  provide  a  fairly  good  approximation  for 
the  profiles  of  the  turbulence  statistics.  On  the  positive  side,  the  analytical  model  described 
herein  affords  a  relatively  simple  method  for  the  prediction  of  canopy  flows  (both  mean  wind 
and  turbulence  statistics)  and  appears  to  capture  enough  of  the  salient  qualitative  features 
of  these  flows  to  be  useful.  With  the  testing  of  the  analytical  model  against  the  urban 
canopy  flows  and  its  modest  success,  it  is  suggested  that  this  analytical  model  for  canopy 
flows  can  probably  be  applied  unmodified  to  a  wide  range  of  canopies,  both  vegetative  and 
urban.  Experience  seems  to  indicate  that  this  imperfect  model  may  give  adequate  solutions 
for  certain  applications. 
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The  one-dimensional  analytical  model  for  canopy  flows  considered  in  this  report  provides  all 
the  statistics  of  the  wind  field  (mean  wind  speed,  kinematic  shear  stress,  velocity  variances, 
and  viscous  dissipation  rate)  required  as  input  to  drive  a  physically-based  model  for  the 
prediction  of  urban  dispersion.  These  flow  statistics  can  be  calculated  rapidly  using  the 
analytical  model  in  comparison  to  the  more  computationally  intensive  second-order  closure 
models  for  canopy  flows.  Taken  as  a  whole,  this  simple  canopy  flow  model  can  probably  be 
used  to  “drive”  a  Lagrangian  stochastic  dispersion  model  for  the  prediction  of  dispersion 
of  materials  released  within  a  canopy.  The  next  logical  step  in  this  work  would  be  to 
construct  a  Lagrangian  stochatic  model  satisfying  the  well-mixed  criterion  that  uses  the 
simple  analytical  model  to  provide  the  input  mean  flow  and  turbulence  quantities  required. 
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Figure  2.  Comparison  of  the  observed  profile  of  the  kinematic  shear  stress 
(open  squares)  in  the  Tombstone  canopy  with  model  predictions. 
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Figure  3.  Comparison  of  the  observed  profile  of  the  turbulence  kinetic 
energy  (open  squares)  in  the  Tombstone  canopy  with  model  prediction. 
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Figure  4.  Comparison  of  the  observed  profile  of  the  streamwise  turbulent 
velocity  standard  deviation  (open  squares)  in  the  Tombstone  canopy  with 
model  prediction  obtained  by  equilibrium  partitioning  of  the  turbulence 
kinetic  energy. 
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Figure  5.  Comparison  of  the  observed  profile  of  the  vertical  turbulent 
velocity  standard  deviation  (open  squares)  in  the  Tombstone  canopy  with 
model  prediction  obtained  by  equilibrium  partitioning  of  the  turbulence 
kinetic  energy. 
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Figure  6.  Comparison  of  the  observed  profile  of  the  streamwise  turbulent 
velocity  standard  deviation  (open  squares)  in  the  Tombstone  canopy  with 
model  prediction  obtained  by  algebraic  stress  model  partitioning  of  the 
turbulence  kinetic  energy. 
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Figure  7.  Comparison  of  the  observed  profile  of  the  vertical  turbulent 
velocity  standard  deviation  (open  squares)  in  the  Tombstone  canopy  with 
model  prediction  obtained  by  algebraic  stress  model  partitioning  of  the 
turbulence  kinetic  energy. 
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Figure  8.  Comparison  of  the  observed  values  of  the  attenuation  coefficient 
at  various  plan  area  indices  for  staggered  and  square  arrays  of  cubical 
obstacles  with  two  proposed  parameterizations. 
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Figure  9.  The  observed  variation  of  the  skin  friction  at  canopy  height  as  a 
function  of  the  plan  area  index  for  staggered  and  square  arrays  of  cubical 
obstacles. 
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Figure  11.  Model  predictions  of  profiles  of  streamwise  turbulent  velocity 
standard  deviation  for  a  square  array  of  cubical  obstacles  with  various 
plan  area  indices.  The  modeled  profiles  were  obtained  using  equilibrium 
partitioning  of  the  turbulence  kinetic  energy.  The  horizontal  bars  show  the 
range  in  the  observations  by  Macdonald  et  al.  [21]. 
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Figure  12.  Model  predictions  of  profiles  of  vertical  turbulent  velocity 
standard  deviation  for  a  square  array  of  cubical  obstacles  with  various 
plan  area  indices.  The  modeled  profiles  were  obtained  using  equilibrium 
partitioning  of  the  turbulence  kinetic  energy. 
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Figure  13.  Model  predictions  of  profiles  of  streamwise  turbulent  velocity 
standard  deviation  for  a  square  array  of  cubical  obstacles  with  various 
plan  area  indices.  The  modeled  profiles  were  obtained  using  algebraic 
stress  model  partitioning  of  the  turbulence  kinetic  energy.  The  horizontal 
bars  show  the  range  in  the  observations  by  Macdonald  et  al.  [21]. 
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Figure  14.  Model  predictions  of  profiles  of  vertical  turbulent  velocity 
standard  deviation  for  a  square  array  of  cubical  obstacles  with  various 
plan  area  indices.  The  modeled  profiles  were  obtained  using  algebraic 
stress  model  partitioning  of  the  turbulence  kinetic  energy. 
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Figure  15.  Model  predictions  of  profiles  of  streamwise  turbulent  velocity 
standard  deviation  for  a  square  array  of  cubical  obstacles  at  high  plan 
area  indices.  The  modeled  profiles  were  obtained  using  equilibrium 
partitioning  of  the  turbulence  kinetic  energy.  Open  symbols  are 
measurements  made  by  Macdonald  et  al.  [21]. 
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Figure  16.  Model  predictions  of  profiles  of  streamwise  turbulent  velocity 
standard  deviation  for  a  square  array  of  cubical  obstacles  at  high  plan 
area  indices.  The  modeled  profiles  were  obtained  using  algebraic  stress 
model  partitioning  of  the  turbulence  kinetic  energy.  Open  symbols  are 
measurements  made  by  Macdonald  et  at.  [21]. 
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